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Abstract. The INTEGRAL satellite extensively observed the black hole binary Cygnus X-l from 2002 November 
to 2004 November during calibration, open time and core program (Galactic Plane Scan) observations. These 
data provide evidence for significant spectral variations over the period. In the framework of the accreting black 
hole phenomenology, the source was most of the time in the Hard State and occasionally switched to the so- 
. called "Intermediate State". Using the results of the analysis performed on these data, we present and compare 

' the spectral properties of the source over the whole energy range (5 keV-1 MeV) covered by the high-energy 

, instruments on board INTEGRAL, in both observed spectral states. Fe line and reflection component evolution 

occurs with spectral changes in the hard and soft components. The observed behaviour of Cygnus X-l is consistent 
with the general picture of galactic black holes. Our results give clues to the physical changes that took place 
in the system (disc and corona) at almost constant luminosity during the spectral transitions and provide new 
measures of the spectral model parameters. In particular, during the Intermediate State of 2003 June, we observe 
in the Cygnus X-l data a high-energy tail at several hundred keV in excess of the thermal Comptonization model 
which suggests the presence of an additional non-thermal component. 

Key words, black hole physics - stars: individual: Cygnus X-l - gamma rays: observations - X-rays: binaries - 
X-rays: general 
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' 1. Introduction soft to contribute significantly. The Thermal Dominant 

State (TDS) instead shows a strong thermal component 

Galactic Black Holes (BH) X-ray binary systems display with a characteristic temperature of at most a few keV 

high-energy emission characterized by spectral and that dom i nates the X-ray spectrum. No, or very weak 

flux variabilities on time scales ranging from millisec- and S p ec t r ally steep, hard X-ray emission is observed; the 

onds to months. These systems are generally found radio emission is quen ched or very faint. This spectrum 

in two major states mainly defined by the relative is generally associa ted with a geometrically thin an d 

variable contributions of soft and hard X-ray compo- optically thick accre tion disc llSha'kura fc Sunvaev 197fl . 

nents, radio spectral properties and timing behaviour In addition to these tw0 canon ical states, other states 

j McClintock ft fiemillard 200.^ | Nowak 2002|) . In the haye been identinedj characterized either by an even 

Hard State (HS), the X-ray and 7 -ray spectrum is greater i uminosity than in the TDS (the "Steep Power 

generally described by a power law model with an i aw State") or by variability and X-ray spectral properties 

exponential cutoff at a few hundred keV, accompanied tly intermediat e between the HS and the TDS 

by relevant radio emission; it can be modeled by thermal fa fi11nni ,,, ; ,|, 199f fr t he "Intermediate State" (IS). 

Comptonization of cool seed photons in a hot electron Cygnug X -1/HDE 226868 is one of the first X-ray 

plasma j Gierlihski et al. lflflfl l Dove et al. 199^1) - The soft binaries detected; it belongs to the BH binary category. 

(~ 0.1-2 keV) black body component is very weak or too Among the brigntest X -ray sources of the sky, it is 

also very variable on different time scales. The assump- 
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tion that Cygnus X-l ranks among the microquasars 
has been confirmed by the detection of a relativistic 
iet JStirling et al. 20011) . Since its discovery in 1964 
(Bo wver et al. 19651) . it has been extensively observed as 
the prototype of BH candidates in radio/optical wave- 
lengths and with all high-energy instruments, from sof t 
X-rays to 7-rays, e.g., with ASCA ijOierliriski et al. "l9 99L 
SIGMA (Salo t ti et al. 1992 ), RXTE 
Dove et al. 1993 IPottschmidt et al. 2003al). BeppoSAX 
Frontera et al. 200lt iDi Salvo et al. 200ll) and CGRO 
McConnell et al. 2000HMcConneil et al. 20021) . This per- 
sistent source, located at ~ 2.4 ± 0.5 kpc (McClintock 
& Remillard 2003, Table 4.1), accretes via a strong stellar 
wind coming from its companion, a 09.71 star of 20 M 
l|Zi61kowski 20051) with an orbital period of 5.6 days. 
The mass function constrains the inclination angle o f the 
system between 25° and 67° ijOierliriski et al. 19991) and 
we adopted the value of 45°. 

Cygnus X-l spends most of its time (90% until 
1998 see, e.g., Gierlinski et al. 1999) in the HS, with 
a relatively low flux in soft X-rays (~ 1 keV) and a 
high flux in hard X-rays (~ 100 keV). Its spectrum is 
roughly described by a power law with a photon index 
T between 1.4-2.1; a break at energies higher than 
~ 50 keV is present. This state is also characterized 
by a large timing variability. Occasionally, the source 
switches to the TDS with V > 2.3. During 1996 June, 
in addition to the dominant black body component and 
the hard component, a h igh-energy tail extendi ng up to 
10 MeV was discovered JMcConnell et al. 2002|) . In this 
state little variability is observed. The IS, in which the 
source exhibits a relatively soft hard X-ray spectrum 
(r ~ 2.1-2 .3) and a moderately stron g soft thermal 
component l|Mendez fc van der Klis 1997j) . often appears 
when the source is about to switch from one main state 
to another. When not associated with a clear spectral 
transition, this state is called a "Failed State Transition" 
(FST). In the IS, the source displays a complicated 
pattern of timing properties. In the past few years, the 
source has been deeply obs erved in the IS and in the 
TPS J/dziarski et al. 20021 fPottschmidt et al. 2003al 
lOleissner et al. 2004al IZHziaxski fc Gierlinski 20041). In 
addition to the thermal and power law components, 
other spectral features can be present in the spectrum: a 
reflection component peaking around 30 keV and, most 
noticeably, a Fe Ka line and Fe edge between 6 and 7 
keV. These features can be visible in both spectral states. 

The INTErnational Gamm a-Ray Astrophysics 
Laboratory (INTEGRAL) mission l|Winkler et al. 20031) 
is an European Space Agency satellite launched on 2002 
October 17, carrying two main 7 -rav instruments, I BIS 
JUbertini et al. 2fwA and SPI JVedrenne et al. 20031) . 
and two X-ray monitors JEM-X llLund et al. 20031) 
Composed o f two detectors, ISGR I JLebrun et al. 20031) 
and PICsIT l|Di Cocco et al. 20031) . the IBIS coded mask 
instrument covers the energy range between 15 keV and 
10 MeV. The SPI telescope works in the 20 keV-8 MeV 
range and the JEM-X monitors provide spectra and im- 
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Fig. 1. RXTE/ ASM daily average (1.5-12 keV) light curve 
of Cygnus X-l from 2002 November to 2004 November 
(MJD = JD - 2400000.5) with the periods of our 
INTEGRAL observations (see text and Tables 1 and 2 
for epoch definitions). 

ages in the nominal 3-35 keV band. As a bright hard X-ray 
source, Cygnus X-l is a prime target for INTEGRAL and 
has been extensively observed during the Performance 
Verification (P V) Phase of the miss i on, when the source 
was in the HS feazzano et al. 20031 iBouchet et al. 200.4 
ICadolle Bel et al. 20041 IPottschmidt et al. 2003hl). 

Pottschmidt et al. (2005) also reported on preliminary 
analysis of the high time-resolution Galactic Plane Scan 
(GPS) observations of Cygnus X-l (up to 2004 April) 
in the 15-70 keV range. In the present work, we report 
the results over a wide energy band (from 5 keV up to 
1 MeV) of several sets of observations of Cygnus X-l, 
including part of the PV-Phase observations not yet 
exploited, the first observations of Cygnus X-l in the 
Open Time program, a larger amount of GPS data 
than previously analyzed and the data from the 2004 
November calibration period. For the first time, up to 1.5 
Ms of INTEGRAL data of Cygnus X-l, collected over 
two years from 2002 November to 2004 November, are 
presented, fully exploiting the broad-band capability of 
all high-energy instruments of the mission. 

2. Observations and data reduction 

Table reports the general periods (epochs) of the ob- 
servations used, giving for each of them the instrument 
data available, date range, exposure (per instrument) and 
observing modes. Epoch 1 includes part of the PV-Phase 
observations of Cygnus X-l. To discuss the time evolution 
of the source, we report here the IBIS/ISGRI light curves 
and hardness ratios obtained during most of the PV-Phase 
observations of Cygnus X-l (first line of Table QJ, but 
since spectral results were presented in previous works, 
we studied more specifically the broad-band spectrum (us- 
ing JEM-X, IBIS and SPI data) only for those PV-Phase 
observations not yet fully exploited, i.e., those performed 
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between 2002 December 9-11 (epoch 1). The Open Time 
observation was performed on 2003 June 7 1 1 (epoch 2) 
with a 5 x 5 dither pattern ll.Tensen et al. 200.51) . The effec- 
tive exposure time was 275 ks for JEM-X2, 292 ks for IBIS 
and 296 ks for SPI. For this latter period, preliminary re- 
sults can be found in Malzac et al. (2004) and Cadolle 
Bel et al. (2004) but we report here the complete study 
of the average spectrum while, in a future work, Malzac 
et al. (submitted) will present the variability properties of 
the source. Epoch 3 and epoch 4 refer respectively to the 
set of Cygnus X-l observations during the core program 
GPS and the 2004 November calibrations. Unfortunately, 
for the short interrupted GPS exposures, the SPI data 
are unusable because a sensitive evaluation of the back- 
ground is not possible: only JEM-X and IBIS/ISGRI data 
are used. The periods of our different INTEGRAL obser- 
vations, presented in Tabled are also indicated in Fig. 
(epochs from numbers 1 to 4). 

We reduced the IBIS and JEM-X data with the 
standard analysis procedures of the Off-line Scientific 
Analysis 0SA 4 . 2 released by the ISDC, whose algorithms 
are described in Goldwurm et al. (2003) and Westergaard 
et al. (2003) for IBIS and JEM-X respectively. A basic 
selection was performed to exclude those pointings too 
close to radiation belt entry or exit, or spoilt by too much 
noise. To avoid uncertainties in the response files for high 
off-axis angles, we also selected the IBIS data of the obser- 
vations for which the source was in the fully coded field of 
view, i.e., with an offset from the telescope axis no larger 
than 4.5°, and JEM-X data for maximum offset angles 
of 3°. Following recommendations of the 0SA 4.2 release, 
IBIS/ISGRI events were selected to have corrected ener- 
gies > 20 keV and rise time channels between 7 and 80. For 
the background correction, we used a set of IBIS/ISGRI 
maps derived in 256 energy channels from empty field ob- 
servations (these maps will be the default IBIS/ISGRI cor- 
rection maps for the 0SA 5 . release) combined with our 
chosen energy bins while, for the off-axis correction maps 
and the response matrices, we used those of the official 
0SA 4.2 release. In the analysis, we considered the pres- 
ence of the two other sources of the region, Cygnus X-3 
and EXOSAT 2030+375, when they were active. For the 
IBIS/ISGRI spectral extraction however, we implemented 
the most recent module (prepared for the 0SA 5 . deliv- 
ery) which is based on the least squares fit done on back- 
ground and efficiency corrected data, using coded source 
zones only. This option minimizes spurious features in 
the extracted spectra, which appear in particular when 
the sources are weak, partially coded and the background 
poorly corrected (A. Gros, private communication). For 
the PICsIT spectral extraction, we took the flux and error 
values in the mosaic image at the best-fit position found 
for the source. We used the response matrices officially re- 
leased with 0SA 4.2, rebinned to the 8 energy channels of 
the imaging output. 

The SPI data were preprocessed with 0SA 4 . 2 us- 
ing the standard energy calibration gain coefficients per 
orbit and excluding bad quality pointings which have 
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Fig. 2. The 20-200 keV IBIS/ISGRI light curve of 
Cygnus X-l from 2002 November 25 until 2004 November 
22 and corresponding HR between the 40-100 and the 
20-30 keV energy bands (see text and Tables 1 and 2 for 
epoch definitions). 

anomalous exposure and dead time values (or with a 
high final y 2 during ima ging). The spiros 9.2 release 
I Skinner & Connell 20031) was used to extract the spec- 
tra of Cygnus X-l, Cygnus X-3 and EXO 2030+375, with 
a background model proportional to the saturating event 
count rates in the Ge detectors. Concerning the instru- 
mental response, version 15 of the IRF (Image Response 
Files) and version 2 of the RMF (Redistribution Matrix 
Files) were used for epoch 1 and 2, e.g., prior to detector 
losses, while versions 17 and 4 respectively were taken for 
epoch 4, e.g., after the failure of two detectors. 

3. Results of the analysis 

As shown in Fig. Q] during the epoch 2 INTEGRAL ob- 
servations, the 1.5-12 keV ASM average count rate of 
Cygnus X-l (~ 1.3 Crab) was larger than during epoch 1 
(~ 290 mCrab) by a factor of 4.5. The derived IBIS/ISGRI 
20-200 keV light curves and Hardness Ratio (HR) of 
Cygnus X-l are shown respectively in Fig. (general 
overview of part of PV-Phase and epochs 1 to 4) and 
Fig. 03 (specific zooms on part of PV-Phase, epochs 1 and 
2). From epoch 1 to epoch 2, while the ASM average count 
rate increased, the 20-200 keV IBIS/ISGRI one decreased 
from ~ 910 to ~ 670 mCrab as shown in Fig. 03 (where, in 
the 20-200 keV range, 1 Crab = 205 cts s" 1 ). This prob- 
ably indicates a state transition between epochs 1 and 2, 
as also suggested by the decrease in the IBIS HR (Fig.|SJ). 
Similar transition, with a change in the ASM light curves 
and an evolving IBIS HR, occurred again during GPS 
data (epoch 3). Figure El (bottom) shows the IBIS HR 
behaviour over the whole 2002-2004 period indicated in 
Table HI 

We sampled epoch 3 in five distinct sub-groups (noted 
a to e) of close pointings which appear to occur, accord- 
ing to Fig. ri] and in different regimes of ASM count 
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Table 1. Log of the Cygnus X-l observations analyzed in this paper. 



Epoch 


Instrument 


Observation Period 


Exposure 


Observation 






(date yy/mm/dd) 


(ks) 


Type 




IBIS 


02/11/25-02/12/15 


810 


staring, 5x5 a , hex " 


1 


IBIS/SPI/JEM-X 


02/12/09-02/12/11 


365/365/31 


5x5 


2 


IBIS/SPI/JEM-X 


03/06/07-03/06/11 


292/296/275 


5x5 


3 


IBIS/JEM-X 


03/03/24-04/09/10 


269/35 


GPS c 


4 


IBIS/SPI/JEM-X 


04/11/22 


8/8/6 


calibration 



Notes: 

a) 5x5 dither pattern around the nominal target location. 

b) hexagonal pattern around the nominal target location. 

c) individual exposures separated by 6° along the scan path, shifted by 27.5° in galactic longitude. 
(Observations indicated on the first line also used, together with epochs 1-4, for Fig. Eland 0) 
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Fig. 3. Top: Zoom on the 20-200 keV IBIS/ISGRI 
light curve of Cygnus X-l from 2002 November 25 until 
December 15 and corresponding HR between the 40-100 
and the 20-30 keV energy bands (average level denoted by 
dashed line). Bottom: Same as above for epoch 2. 



rate and of average IBIS HR. The data of each epoch 
(and sub-group) were summed to obtain an average spec- 
trum studied separately. We added 3% systematic er- 
rors for JEM-X (5-30 keV range), IBIS (20 keV-1 MeV 
range) and SPI (22 keV-1 MeV range) and fitted the 



resultant spectra si multaneously using XSPEC vl 1 . 3 . 
l)Arnaud et al. 199^ . In order to account for uncertain- 
ties in the cross-calibration of each instrument, a multi- 
plicative constant was added in the spectral fits to each 
instrument data set: it was set free for IBIS and SPI and 
frozen to 1 for JEM-X. 

Several models were used in the course of the spec- 
tral analysis. In XSPEC not ation, we used a mu lticolour 
disc black body diskbb iMitsuda et al. 19841) p lus a 
Comptonization model COMPTT |Titarchuk 19941) and, 
when necessary, we added a Gaussian for the Fe line 
with the GAUS SIAN model and the reflect ion compo- 
nent REFLECT llMagdziarz & Zdziarski 19951) . This lat- 
ter component models the X-ray reflection of the comp- 
tonized radiation from neutral or partially ionized mat- 
ter, presumably the optically thick accretion disc (Done 
et al. 1992, Gierlihski et al. 1997, 1999, Magdziarz & 
Zdziarski, 1995). For this source, we always used a 
fixed absorption column densi ty Nn of 6 x 10 21 cm~ 2 
l|Ba,luciriska,-Church et al. 199.51) . We also tied the input 
soft photon temperature kTo of the COMPTT model to the 
inner disc temperature kT- m value found by the diskbb 
model. In order to compare all our data (from epochs 1 
to 4) with the same model, we show the parameters ob- 
tained from the current fitted model described above (mul- 
ticolour disc black body plus Comptonization convolved 
by reflection and Gaussian when necessary) . We also tried 
more complex models such as COMPPS and EQPAIR, de- 
veloped respectively by Poutanen & Svensson (1996) and 
Coppi (1999), coupled to the GAUSSIAN and the DISKBB 
(or diskpn for eqpair only, see Section 3.2) models: we 
present such results only for epoch 2, when the statistics 
were significantly better, the instrument configurations 
stable and during which the presence of a non-thermal 
component appeared more pronounced. 

3.1. The hard state spectrum 

Figure Q] shows the resultant count spectrum obtained 
in epoch 1 (2002 December 9-11) from 5 keV up to 
1 MeV with JEM-X, SPI, IBIS/ISGRI and PICsIT data. 
A simple power law model does not fit the spectra well 
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Table 2. Best-fit parameters of Cygnus X-l for the current thermal model in the different observation epochs. 



Epoch Dates Disc Norm. a fcT in or kT kT c r E Fc line n/2n" xLd 
(MJD) (keV) (keV) (keV) (dof) 

1 52617-52620 - 0.20 (frozen) 67 t % 1.98 t u a ^ - 0.25 1.45 (230) 

2 52797-52801 250 ±|| 1.16 ± 0.07 100 tf 7 0.98 tlH 7.07 lg;i? 0.57 ±g:g§ 1.69 (236) 
3a 52710-52780 - 0.20 (frozen) 68 ±f 2 2.08 t° f 4 6.48 ± 0.13 0.32 ig;J$ 1.07 (190) 
36 52801-52825 312 1-15 ± 0.03 93 ± 42 0.80 1^ 6.40 ± 0.73 0.58 ±o.is °- 93 ( 19 °) 
3c 52990 361 ±f 7 0.99 ± 0.08 58 ±\\ 1.60 Igjo 6.96 ± 0.19 0.23 ±%j& 0.99 (190) 
3d 53101-53165 0.20 (frozen) 56 ± 12 2.28 ±om 6 - n ± °- 2 6 0.27 ± 0.06 0.81 (190) 
3e 53240-53260 132 ± 10 1.39 ± 0.77 48 ± 2 g 1.85 ±° .m 6 - 49 ± 0.38 0.49 1.56 (190) 
4 53335 232 tf 2 Ijj (frozen) 128 ±63 0-74 ±g;|| 7.78 ±g^ 0.47 ±g^g 0.97 (221) 



Notes: 

a) Disc normalization K is given by K = (R/D) 2 cos 9 where R is the inner disc radius in units of km, D is the distance to 
the source in units of 10 kpc and 9 the inclination angle of the disc. 

b) Solid angle of the reflection component. 

Model applied in XSPEC notations: constant*wabs*(diskbb+gaussian+reflect*comptt) with Nn fixed to 6 x 10 21 cm~ 2 
and kT value tied to disc kT in . Errors are at 90% confidence level (A^ 2 = 2.7). 
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Fig. 4. Spectra of Cygnus X-l in 2002 December 9-11 
(epoch 1) with the JEM-X (black), SPI (green) and IBIS 
(ISGRI: red; PICsIT: blue) data. The best-fit model is 
a Comptonization model with reflection (see Table 2). 
Residuals in a units are also shown. 



(photon index T of 1.9 ± 0.1 and reduced chi-square 
xfed = 12-90 with 213 degrees of freedom, hereafter 
dof). A cutoff in the model, with a typical folding en- 
ergy of approximately 150 keV, clearly improves the fit 
{X 2 red = 2.12 with 212 dof) and better describes the 
available data. Since a cutoff power law is usually at- 
tributed to thermal Comptonization, we replaced this 
phenomenological model by a more physical model of 
Comptonization (COMPTT, Titarchuk 1994). Some resid- 
uals were still visible around 10 keV so we added a model 
of reflection (with an inclination angle equal to 45°) to 
account for this excess. The final tested model there- 
fore includes thermal Comptonization convolved by re- 
flection (reflect), with solar abundances for Fe and He 
l|Anders fc Grevesse 198fl(l . We obtain a plasma tempera- 
ture kT e of 67 keV with an optical depth r of 1.98 and 
fl/2n = 0.25, with xled = L45 (230 dof). The disc black 
body is very weak or below the energy range of JEM- 
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Fig. 5. Epoch 1 unabsorbed EF(E) spectrum of 
Cygnus X-l along with the best-fit model described 
in Table 2 with the JEM-X, SPI and IBIS (ISGRI 
and PICsIT) data. Dotted: Reflection. Long Dashes: 
Comptonization. Thick: Total model. 

X: this component was not used in our fits. As it gives 
no contribution, we froze the kTo temperature of COMPTT 
at 0.20 keV. Normalization constants between instruments 
(JEM-X, IBIS/ISGRI, SPI, IBIS/PICsIT) are respectively 
equal to 1, 1.2, 1.3 and 0.9. 

Note that some residual points are visible in this plot, 
in particular for the SPI data at low energies (and for 
the IBIS data around 200 keV). This is mostly related to 
the non perfect cross-calibration between the INTEGRAL 
instruments. Indeed, by fitting separately the different 
instrument data, the residuals are reduced, with little 
change in the spectral parameters. Future improvement in 
the cross-calibration of INTEGRAL telescopes will allow 
a better determination in the relative flux normalizations 
and also a better agreement of the derived spectral shapes. 
The IBIS configuration was not stable in the first phases of 
the mission: in particular, the PV-Phase spectra may suf- 
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Fig. 6. Spectra of Cygnus X-l in 2003 June 7-11 
(epoch 2) with the JEM-X (black), SPI (green) and IBIS 
(ISGRI: red; PICsIT: blue) data. The best-fit model is a 
multicolour disc black body and a Comptonization model 
with Gaussian and reflection components (see Table 2). 
Residuals in tr units are also shown. 



fer from the fact that the IBIS responses were built from 
the Crab nebula observations performed in later periods, 
with a defined and stable configuration. 

In spite of these caveats, all spectra overlap relatively 
well and do define the same set of spectral parameters. 
Figure El shows the resultant EF(E) spectrum and its 
best-fit with the JEM-X, IBIS and SPI data. Table 
summarizes all fit results. While the 20-100 keV lumi- 
nosity is 6.5 x 10 36 ergs s _1 (with a distance to the 
source fixed at 2.4 kpc), the bolometric luminosity (ex- 
trapolated from 0.01 keV to 10 MeV) has the value of 
2.2 x 10 37 ergs s _1 . The best-fit parameters we ob- 
tain are consistent with those found i n BH binaries in 
the HS JMcClintock fc Remillard 20031) and with previ- 
ous results reported on Cygnus X- l in the HS, in par- 
ticular from the INTEGRAL data JBazzano et al. 200.4 
iBouchet et al. 2003UPottschmidt et al. 2003hl). These au- 
thors generally give plasma temperatures between 50-100 
keV and optical depths in the 1.0-2 range. In particular, 
our SPI spectrum is fully compatible with the one reported 
by Bouchet et al. (2003). These authors observed that the 
derived spectrum had an excess component with respect 
to a Sunyaev & Titarchuk (S-T) Comptonization model 
(1980) and that it was better described by a cutoff power 
law with r = 1.5 and a cutoff energy of 155 keV. We find 
exactly the same results for our SPI spectrum alone if 
we fit it with an S-T model or a cutoff power law. This 
excess emission, with respect to the simple S-T model, 
was already observed in a number of cases for the HS 
spectra of BH systems, e.g., Jourdain & Roques (1994). 
However, with the Titarchuk (1994) model (which better 
describes the Comptonization for high plasma tempera- 
tures), no significant excess is observed in epoch 1 spectra 
of Cygnus X-l. 



3.2. The transition to a softer state in 2003 June 

As visible in the change of the ASM light curve of 
Cygnus X-l (see Fig. and its corresponding high- 
energy HR shown in Fig. (average value of 1.1 for epoch 
2 compared to more than 1.4 for epoch 1), the source 
softens during epoch 2. Although a single power law does 
not fit the data properly {xl e d = 23.51 with 245 dof), 
the derived slope is softer (r = 2.2 ± 0.1). As for the HS 
spectrum analysis, an inspection of the residuals reveals 
what kind of component can be added to build a proper 
model. The need for each new component is then checked 
with the results of the fits and the final model is further 
verified by repeating the procedure with components 
added in a different order. 

First, we included a cutoff in the power law and 
we improved the x 2 by A% 2 = 7.8: this indicates that 
the new component is significant at more than the 95% 
confidence level. A cutoff in the power law model, with a 
typical folding energy of approximately 200 keV, better 
describes the available data, so we tried the physical 
Comptonization model (comptt) as we did for epoch 1. 
Moreover, since very large residuals are visible in the soft 
X-rays and since we suspect a transition to the IS, we 
added a multicolour disc black body (diskbb), which is 
in fact required by the data (xled = 14-87 and 239 dof 
without this component). Some residuals around 10 keV 
indicate the need for a reflection component. Therefore, 
the best model for the continuum consists of a multicolour 
disc black body (diskbb) and thermal Comptonization 
(comptt) convolved by reflection (REFLECT, with pa- 
rameters as above). The Comptonization temperature 
kTo was fixed to the kT- m value returned by the diskbb. 
Finally, we added a Gaussian to account for the residuals 
in the JEM-X data around 6-7 keV (xled = 2 - 42 and 237 
dof without this component). 

Table summarizes the best-fit parameters and the 
xl e d obtained from 5 keV up to 1 MeV. We obtained a 
plasma temperature kT c of 100 1 17 keV and an optical 
depth t of 0.98 tp^g, respectively higher and lower than 
in epoch 1 (kT e = 67 keV, r = 1.98 toil)- The 
inner disc temperature reached 1.16 keV and a significant 
line is detected at a centroid energy of 7.07 keV, with 
an Equivalent Width (EW) of 1.4 keV. With the same 
assumed distance, the luminosity is 6.5 x 10 36 ergs s _1 
in the 0.5-10 keV range and 5.2 x 10 36 ergs s -1 in the 
20-100 keV band. The bolometric luminosity, extrap- 
olated from 0.01 keV up to 10 MeV, has the value of 
2.0 x 10 37 ergs s _1 ; the disc accounts for 26 % of the 
total luminosity. 

The derived disc normalization is possibly underes- 
timated and not well constrained by the JEM-X data 
(which start at 5 keV) . The multicolour disc black body is 
only an approximation of the soft component: the direct 
derivation of physical parameters of the disc from the 
best-fit values may suffer from some important effects 
I Merloni et al. 200(1) . According to these authors, the 
dominant effect seems to be that, in the inner part of 
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the disc, the opacity is dominated by electron scattering 
rather than free-free absorption. The net result is that 
the derived temperature given by the kT- ln parameter 
overestimates the effective inner temperature by a factor 
of 1.7 or more l|Shimura fc Takahara, 1995^ . This has 
an effect on the estimation of the inner disc radius as 
well. Assuming an inclination angle of 45°, the derived 
internal radius from the best-fit disc normalization has 
the unphysical value of 4.50 1Lq' 7 5 q km. For a 10 M Q 
BH, this value corresponds to 0.15 R s (where R s is the 
Schwarzschild radius), i.e., smaller than the last stable 
orbit even around a maximally spinning Kerr BH. 

A comparatively strong reflection component 



0.57 



+0.09 
-0.06 



while it was 0.25 



+0.03 

-0.04 



for epoch 1) is also necessary to fit the data 
(Xred = 4-57 and 237 dof without this component). 
Our £V2tt value is co nsistent with ASCA RXTE 
ijGierlinski et al. 199ft IGierlinski et al. 19991) and 
BeppoSAX observations ijFrontera, et al. 2001]) . The 
best-fit model, over-plotted on data from JEM-X, SPI, 
IBIS/ISGRI and PICsIT, is reported in count units in 
Fig. El Figure[3shows the resultant EF(E) spectrum with 
its best-fit. Note that we have rebinned, for illustrative 
purposes, the SPI high-energy points above 750 keV to 
reach the level of 3er, but the fit was performed on the 
original energy channels of the spectrum (50 bins over 
the 22 keV-1 MeV band). Normalization constants are 
respectively equal to 1, 1, 1.2 and 1. Residuals again show 
that improvements in the instrument cross-calibration 
and responses are still needed to obtain a fully satisfac- 
tory fit of the INTEGRAL combined data, but we believe 
that the general spectral model and parameters are well 
determined. 

Considering the behaviour of the ASM light curve 
(FigHJ, the evolution of the high-energy IBIS HR (Fig.0 
lower panel) , the relative softness of the spectrum and the 
presence of a relatively strong hard energy emission, it ap- 
pears that during the 2003 June observations Cygnus X-l 
was in the IS (or in the FST). This is also confirmed 
by radio observations of Malzac et al. (2004) who sug- 
gested that the fluctuations of the radio luminosity were 
associated with a pivoting of the high-energy spectrum. 
The derived thermal Comptonization parameters are 
consistent with those found in BH binaries in soft states 
jMcClintock fc R,emillard~20 03). 

As it can be seen in Fig. an excess with respect to 
the Comptonized spectrum above 400 keV is observed 
in the SPI data that is not present in epoch 1. We tried 
different models of the background, we used different 
selections of pointings with two data processings (the 
standard pipeline and the specific softwares developed 
by the SPI instrument team); however, this feature 
is always present. To evaluate the possible presence 
of an instrumental feature, we re-analyzed the Crab 
observation performed with SPI on 2003 February 22-24, 
with the same configuration and dither pattern as those 
of our Cygnus X-l observations. We obtained a spectrum 
consistent with the one reported by Roques et al. (2003). 



With 3% systematics added to the 200 ks SPI data, 
the Crab spectrum is well described by a power law 
with T = 2.1 over the 22-1000 keV energy range and no 
high-energy excess above 400 keV is observed. This test 
shows that it is unlikely that the observed high-energy 
tail in the Cygnus X-l IS spectrum is due to a systematic 
effect. 

We fitted this excess with an extra component, in 
addition to our current Comptonization model. Using a 
power law for all the data (except JEM-X), we obtain 
a best-fit photon index T of 2.12 ^o'lg. The improve- 
ment in the xled ls relatively important (1.65 with an 
absolute chi-square reduced by 10 for 2 additional free 
parame ters), making this compon ent significant. An 
F-test i Bevington fc Robinson 1992T) provides a chance 
probability that this improvement is due to fluctuations 
of about 10 -2 (but see Protassov et al. 2002 about 
the limitations of the F-test). At the same time, the 
parameters obtained for the current model are different 
to those obtained without the power law. We obtained a 
plasma temperature of 50 keV, an optical depth of 2.1 and 
fl/2ir is around 0.49. While the IBIS high-energy data 
do not require this additional component, IBIS points 
are also consistent, due to the large error bars, with this 
additional power law model. A better Xred ( e Q ua l to 1.55, 
with an absolute chi-square reduced by 35, giving now a 
chance probability of 10 -5 for the additional power law 
that the improvement is due to fluctuations) is obtained 
when PICsIT data are neglected. Since the SPI telescope 
is certainly the best calibrated INTEGRAL instrument 
at energies higher than 400-500 keV, we consider this 
feature to be significant. 

Consequently, in order to account for this high-energy 
emission with more physical models, we fitted the data 
with the hybrid models COMPPS and eqpair (coupled 
to the usual disc and Gaussian line components). They 
combine both thermal and non-thermal particle distribu- 
tions in the calculation of the emergent spectrum, as fully 
described by Poutanen & Svensson (1996) and Coppi 
(1999). Those models use a hot plasma cloud, mainly 
modeled as a spherical corona around the compact object 
(configuration taken in our fits) or as the well-known 
"slab" geometry, illuminated by soft thermal photons com- 
ing from an accretion disc. These photons are Compton 
scattered by both thermal (Maxwellian) and non-thermal 
(power law) electrons that lose energy by Compton, 
Coulomb and bremsstrahlung interactions. The number 
of electrons is determined by r, the corresponding (total) 
vertical Thomson optical depth of the corona. 

Both EQPAIR and COMPPS models allow one to inject 
a non-thermal electron distribution with Lorentz factors 
between 7 m ; n and 7 max and a power law spectral index T p . 
In the eqpair model, the system is characterized by the 
power Li supplied to its different components, expressed 
as the dimensionless compactnesses U = (LiCx) / Rm e c 3 
(where R is the corona size). This model has two advan- 
tages with respect to the COMPPS: first, it is valid both for 
high kT c and low r and vice- versa (while COMPPS should 
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be used for t < 3 only). Secondly, the EQPAIR model 
also takes into account the electron-positron annihilation 
process. We used the diskpn and diksbb routines in 
the XSPEC package to describe the soft emission from 
the accretion disc in the two hybrid models (eqpair and 
then COMPPS). With the DISKPN, the main characteristic 
is that the disc spectrum incident on the plasma is 
computed assuming a pseudo-Newtonian potential for 
the accretion disc, extending from 3 (the minimum stable 
orbit) to 500 R s . 

Both the purely thermal and the thermal/non-thermal 
best-fit parameters of these hybrid models for the Cygnus 
X-l spectrum of 2003 June are reported in Table The 
pure thermal models obviously do not improve the fit 
with respect to the thermal Comptonization model, but 
when the effect of a non-thermal electron distribution is 
included, the x 2 are reduced by a significant amount and 
the high-energy component appears well fitted by these 
models, in particular by the non-thermal eqpair. 

First, for eqpair, with frozen values of 7 m i n and 7 max 
at respectively 1.5 and 1000 and with the same value for 
iVn and abundances as before, we obtain the dimension- 
less parameters lh/l s = 4.57 1q!s7 an d lnth/lh = 0.16 ^o'cm 
(where l s , l t h, Inth and l n = 2th + ^nth correspond to 
the compactness in soft disc photons irradiating the 
plasma, in direct thermal electron heating, in electron 
acceleration and in total power supplied to electrons in 
the plasma respectively). These parameters indicate that 
the power in the non-thermal component is ~ 16% of 
the total power supplied to the electrons in the corona. 
With a Xred = 1-55 (232 dof), clearly better than the 
current epoch 2 model and the pure thermal version of 
the EQPAIR model (values also reported in Table EJ) , the 
derived thermal values of t, fl/2ir, E-p c centroid and EW 
match, within the uncertainties, the parameters obtained 
in Table [2j The value of kT c decreases from the pure 
thermal model as expected and it is similar to the value 
we get when adding a power law to our current model. 
While fitting a soft state spectra from 0.5-200 keV with 
the same hybrid model, Frontera et al. (2001) found T p 
equal to 2.5 ± 0.1, an EW - 300 eV and Q/2ir around 
0.63. These values match well with our non-thermal 
results. However, they obtained a lower lh/U ~ 0.36 and 
a higher l n th/lh ~ 0.77. For what they called a typical 
soft state, McConnell et al. (2002) found a T p compatible 
with ours (around 2.6) but, again, lh/l s ~ 0.17 and 
lnth/lh ~ 0.68 different from our values. For another soft 
state, Gierlihski et al. (1999) found even higher values 
for these latter three parameters, in the range 2.6-3.4, 
0.25-1.6 and 0.77-1. All these comparisons suggest that, 
in epoch 2, Cygnus X-l was not in the typical soft state, 
but was rather in the IS. 

Second, the thermal part of the hybrid model COMPPS 
fits the data for a xled = 1 - 72 ( 233 dof )- Most of the 
crucial parameters like kT e , the optical depth, fl/2ir and 
the -Efc centroid are compatible with those obtained with 
our current model (reported in Table EJ with a compa- 
rable xled- TT^g to improve our fits, we then used the 



non-thermal part of the COMPPS model with an electron 
distribution of index 3.72 io'22 and we re P or t all the 
derived parameters in Table 03 Comparing the different 
xled (now equal to 1.69 for 230 dof), the improvement is 
significant but not as strong as the eqpair non-thermal 
model. Again, kT c decreases from the value obtained 
with a pure thermal model, which is compatible with 
McConnell et al. (2002). Using the non-thermal part of 
the COMPPS model to fit a soft state, these authors found 
temperatures (~ 60 keV), reflection factors and energy 
index T p (~ 3.5) again compatible with ours. 

Some authors, e.g., Belloni et al. (1996) and Ibragimov 
et al. (2005, submitted) claim that the soft states ob- 
served in 1996, 1998 and 1999 were in fact IS. Modeling 
the spectra with eqpair and a non-thermal high-energy 
component, Ibragimov et al. showed spectral results 
consistent with the non-thermal values we report here, 
including a high l n /l s value. According to the different 
soft states of Cygnus X-l analyzed, they obtained for this 
parameter values between 4.2 and 12 and for r, values in 
the range 0.9-1.6. These results are consistent with ours. 

We find that the non-thermal hybrid EQPAIR model 
the best able to fit our epoch 2 data: this model clearly 
accounts for the high-energy tail observed. The non- 
thermal component represents 16% of the total power 
supplied to the electrons in the corona and the inferred 
luminosity in the 20-100 keV range is 6 x 10 36 ergs s _1 
while the bolometric one is 3.3 x 10 37 ergs s _1 . 
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Fig. 7. Epoch 2 unabsorbed EF(E) spectrum of 
Cygnus X-l along with the best-fit model described in 
Table 2 with the JEM-X, SPI and IBIS (ISGRI and 
PICsIT) data. Dashed: Disc. Dotted- Dashed: Gaussian 
line. Dotted: Reflection. Long Dashes: Comptonization. 
Thick: Total model. 
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Table 3. Best-fit parameters of Cygnus X-l for thermal/non-thermal hybrid models in the epoch 2 observations. 



Model 


eqpair (thermal) 


eqpair (non-thermal) 


compps (thermal) 


compps (non-thermal) 


l}x/ls 


5.3 +"■" 

— 0.4 


4.57 + r 87 

— U.esY 






/n-f-Vi / /li 
•'IlLIl / "11 




16 +S'ii 

W-J - U —0.08 






T 


n +0.02 
(J. DO n no 

— o.oy 


0.49 +°f 2 

— 0.0*i 


0.57 +°'°i 

— 0.01 


1.43 ± 0.10 


r„ 




2.4 ±1:1 




o 79 +0.17 
— U./z 






1.5, 1000 (frozen) 




1.34, 1000 (frozen) 


n/27r 


n 7 o +0.07 
u - '° -0.03 


0.63 ± 0.08 


0.60 ± 0.05 


0.87 +_° it 


E Fc (keV) 


7 IS + 016 
'•18 -0.14 


7.06 ± 0.06 


7 20 + 10 


7.14 ± 0.09 


EW (eV) 


287 


379 


238 


412 


kT (keV) 


1 44 +0.07 

-0.03 


1 Sq +0-27 


1.20 (frozen) 


1.20 (frozen) 


kT e (keV) 


103 a 


42 a 


111 ±1 


39 If 


Xred (dof) 


1.68 (234) 


1.55 (232) 


1.72 (233) 


1.69 (230) 



Notes: 

a) The electron temperature is calculated for the best-fit model (i.e., not free parameters). 

Model applied in XSPEC notations: constant*wabs*(diskpn+gaussian+eqpair) or con- 
Stant*wabs*(diskbb+gaussian+compps) with JVh fixed to 6 x 10 21 cm -2 and kTo value tied to disc kTi n . Errors are 
at 90% confidence level (A x 2 = 2.7). 



3.3. Evolutions occurring during Epochs 3 and 4 
(2003-2004) 

We separated epoch 3 in five sub-groups as follows: a) be- 
fore epoch 2 (MJD - 52710-52780); b) just after epoch 

2 (MJD ~ 52801-52825); c) at MJD = 52990 when the 
source switches again to a softer state; d) when a harder 
state is then observed (MJD ~ 53101-53165); e) obser- 
vations when the ASM count rate slightly increases while 
the IBIS HR decreases (MJD - 53240-53260). The IBIS 
HR shown in Fig.[2]indicates spectral evolution with time: 
from groups a to e, the mean HR value is 1.55, 0.8, 1.1, 
1.4 and 1.05. Accordingly, a simple power law model fits 
the data with different indexes: T = 1.9, 2.3, 2.2, 1.9 and 
2.3 ± 0.1 for groups a to e. 

We fitted from 5 keV to 400 keV the JEM-X and 
IBIS/ISGRI spectra of the data collected during these 
sub-groups, using our current Comptonization (comptt) 
model for the HS and the IS defined above. We fixed 
the seed photon temperature kTo to the kT ln value given 
by the diskbb (or frozen to 0.20 keV when the disc 
emission is not detected). The best-fit parameters are 
listed in Table Normalization constants between the 
two instruments are very close to each other (~ 1-1.2 
for IBIS/ISGRI when JEM-X constant is frozen to 1). 
Variations in amplitude and overall spectral shape are 
observed between the epoch 3 sub-groups. While no disc 
emission is detected for groups a and d, a disc component 
is required for groups 6, c and e. The EWs are also chang- 
ing: the softer the source, the larger the EW (~ 110 eV for 
group a compared to ~ 780 eV for group e) . These results 
and the IBIS HR (Fig.0 indicate that, during sub-groups 

3 a and d, Cygnus X-l was in a HS (as in epoch 1) while, 
in sub-groups 3 b, c and e, the source was in a softer state. 

During epoch 4, the simple power law slope is 
T = 2.2 ± 0.1 and the mean IBIS HR value is 1.1. 
This indicates that the source was again in a softer state 
than the HS. The count spectrum is presented in Fig. El 



Figure shows the resultant EF(E) spectrum and its 
best-fit model overplotted on the data obtained from 
JEM-X, IBIS/ISGRI and SPI (high-energy IBIS points 
are rebinned at the 3a level above 200 keV for illustrative 
purposes). Cross-calibration constants are respectively 1, 
1 and 1.3. Using the current Comptonization (comptt) 
model defined above for epoch 2, we present in Table[2]the 
best-fit parameters obtained. Since we obtained a best-fit 
solution similar to the one of epoch 2, we fixed, to de- 
fine the errors for the other parameters, the kT ln to the 
epoch 2 value of 1.16 keV and we obtained kT c = 128 keV, 
t = 0.74, Q/2tt = 0.47 to.u ( close to the e P och 2 value of 
0.57 ±o° l) and E Fc = 7.78 keV. Althought the tempera- 
ture is little constrained and errors are large, these values 
indicate that the source was again in a rather soft state. 
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Fig. 8. Epoch 4 spectra of Cygnus X-l with the JEM- 
X (black), SPI (green) and IBIS/ISGRI (red) data along 
with the best-fit current model (see Table 2). Residuals in 
a units are also shown. 
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Fig. 9. Epoch 4 unabsorbed EF{E) spectrum of 
Cygnus X-l along with the best-fit model described 
in Table 2 with the JEM-X, IBIS/ISGRI and SPI 
data. Dashed: Disc. Dotted- Dashed: Gaussian line. Dotted: 
Reflection. Long Dashes: Comptonization. Thick: Total 
model. 

4. Discussion 

During the broad-band (5 keV-1 MeV) INTEGRAL 
observations of Cygnus X-l presented in this paper, the 
source was detected in at least two different spectral 
states. For epoch 1 (part of PV-Phase) and during some 
of the GPS observations (epoch 3, sub-groups a and 
d), the source was in the typical HS with a high-energy 
spectrum extending up to 800 keV, well characterized by a 
thermal Comptonization model. Parameters are typical of 
BH binaries and consistent with previous observations of 
Cygnus X-l in the HS. We find that the Comptonization 
component must be modified by reflection. The sub- 
tended angle we derived is also compatible with values 
previously found for this source in the HS, either from 
Ginga and CGRO/ OSSE (kT c ~ 100 keV, r ~ 1-2, 
and fl/2iT ~ 0.19-0.34, Gierlihski et al. 1997), or from 
BeppoSAX (fcT c - 60 keV, r ~ 1-2 and Q/2n = 0.25, 
Frontera et al. 2001). Di Salvo et al. (2001) observed even 
higher temperatures for the HS, from 111 to 140 keV, 
and S1/27T ~ 0.1-0.3, within our error bars. Also, our 
results are compatible with those reported by Bazzano 
et al. (2003), Bouchet et al. (2003) and Pottschmidt 
et al. (2003b) who used INTEGRAL data from other 
PV-Phase observations, when Cygnus X-l was also in 
the HS. We did not detect any Fe line nor disc emission 
during this period. Indeed, we could determine a 90% 
confidence level upper limit of 94 eV EW for a narrow 
line (a = 0.1 keV) at 6.7 keV and of 172 eV EW for a 
broad one (a = 1 keV). These limits are compatible with 
the typical EWs (~ 150 eV) of broad lines o bserved in 
HS from this BH system l|Frontera et al. 200lh . The lack 
of strong Fe line and disc emission is not surprising since, 
in the HS, the disc does not extend close to the BH, its 
inner disc temperature is low and contribution at > 5 keV 



is negligible. 

Cygnus X-l underwent a clear evolution to a softer 
state from epoch 1 to epoch 2. We indeed observed 
in 2003 June INTEGRAL data hardness variations 
(Fig. and El and photon index changes, along with the 
appearance of a significant soft component well fitted by 
a disc black body model. The changes in the high-energy 
component are obvious from the best-fit parameters of the 
current model reported in Table The Comptonization 
parameter y is characterized by both values (kT c , r) 
where y = (fcT c /m c c 2 ) max (r,r 2 ). While y ~ 0.51 t£j$ 
for epoch 1, its value drops to 0.19 ^06 f° r epoch 2. 

As shown previously, these changes were correlated 
with evolutions of the RXTE/ ASM light curve (Fig. [TJ 
which indicate the rise of the soft disc emission. The 
increase in disc emission (up to about 26% of the 
bolometric luminosity) combined with an increase of 
the inner disc temperature to about 1.16 keV suggest 
that the accretion disc has extended down very close 
to the BH horizon. However, the inferred inner radius 
(4.5 km) is not compatible with the size of the innermost 
stable circular orbit, even for a maximally spinning 
Kerr BH. Besides uncertainties on the inclination angle 
and on the disc normalization, a number of effects can 
lead to such unphysical estimations. As mentioned in 
Section 3.2, electron scatterings can pr oduce a diluted 
black body ijShiimira, & Takaha.ra 19951) and cause the 
observed temperature kT in to be higher and the disc 
normalization lower. Consequently, the inner disc radii 
could be underesti mated by a factor of 5 or more 
jMerloni et al. 2000l) . Another possible explanation 
for such a high temperature could be that the disc 
emission is Comptonized by a different warm corona 
at a temperature of a few keV, as some times suggested 
( Frontera et al. 200lilMalzac et al. 20051) . 

The softening of the hard component and the appear- 
ance of the disc emission were accompanied by significant 
increase in reflection and changes in Fe line emission. 
During epoch 2 data, the Fe line is needed in the spectra 
(while not for epoch 1) and fl/2iT increases (it is more 
than twice the value found for epoch 1): this again 
indicates that a larger disc, closer to the BH, reflects 
more radiation coming from the hot plasma. Figures 03 
and show the resultant unabsorbed EF(E) spectra 
for both epochs 1 and 2: we clearly see differences in 
amplitudes and in the relative contributions of the various 
components, from soft to hard X-rays/7-rays. Similar 
spectral transitions are observed during epochs 3 and 4: 
y = 0.57 for epoch 3 sub-group a (while it drops to 
0.14 ±o Jg for group 6), and y = 0.18 for epoch 4. 

The y evolution, combined with the changes in the IBIS 
HR and in the ASM light curve, indeed shows transitions 
from the HS (epoch 1, epochs 3 a and d) to softer states 
(epochs 2, 3 6, c, e and 4). 

The transition to a soft state is generally attributed 
to an increase in the accretion rate. In this case, the total 
luminosity should increase significantly. For example, the 
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total luminosity measured by Frontera et al. (2001) in the 
TDS of Cygnus X-l was about a factor of 3 higher than 
in the HS, with a contribution of the disc emission higher 
than 50% of the total luminosity. From our analysis, 
the spectral transition seems to occur with only a slight 
change in the bolometric luminosity, mainly due to the 
appearance of a high-energy tail which represents less 
than 2% of the Eddington luminosity (~ 1.5 x 10 39 ergs 
s _1 for a 10 M Q BH). Therefore, we can interpret the 
component evolution as a spectral pivoting between soft 
and hard X-rays at almost constant luminosity. This 
has already been reported for ot her BH such as for th e 
microquasar XTE J1550-564 ^Rodriguez et al. 20031) . 
When the total luminosity does not change much, the 
spectral evolution cannot be explained simply by a large 
variation in the accretion rate. An additional physical 
parameter (for example linked to the temperature and 
geometry of the Comptonizing cloud or to the magnetic 
field) must probably vary in the system to trigger the 
spectral changes from the HS to a softer state at almost 
constant luminosity. 

While our data start at 5 keV, leading to a possible 
underestimation of the bolometric luminosity (which in a 
real soft state is dominated by the soft component), the 
near constancy of the luminosity between epochs 1 and 
2 and other facts show that the soft state we observed 
in this latter period is not the typical TDS observed in 
Cygnus X-l. Our spectral analysis and the comparisons 
with previous results obtained on Cygnus X-l (see 
Section 3.2) led us to conclude that, during epoch 2, the 
source was rather in the IS. The need for a cutoff in the 
high-energy spectrum of epoch 2 shows that the thermal 
Comptonization is still a dominant process for the 
high-energy component. The hybrid thermal/non-thermal 
models lead to a significant improvement in the fit of our 
data, compared to the thermal Comptonization models, 
like COMPTT, or the pure thermal versions of eqpair and 
COMPPS. However, our best-fit non-thermal parameters 
are sometimes intermediate between hard and soft state 
values (like lh/l s , ^nth/^h or T p ): this again shows that 
Cygnus X-l was in the IS (or in the FST). Such a 
conclusion was also reached by Malzac et al. (2004) 
on the basis of the observed correlation of the spectral 
hardness with the radio flux during this period. 

In recent years, it has become apparent that in 
the HS, the BH binaries become bright in radio and 
display clear correlations betwee n the X-ray and r adio 
luminosities l|Corhel et al. 20031: iGallo et al. 20031 ) as 



often observed for Cygnus X-l jBrocksopn et al. 1 99f. , 
iGleissner et al. 2004ri iNowak 200.51) . Models where the 
base of a compact jet plays a major role in the phys- 
ical processes of such BH systems have been proposed 
(Markoff et al. 2001). In such scenarios, the high-energy 
emission seen during the HS is interpreted as syn- 
chrotron emission from the jet that extends from radio 
to hard X-rays, naturally explaining the correlations 
observed during the HS. Those models explain the 
observed (or inferred) outflows: the radio emission is 



proportional to the jet power which in turn is correlated 
with the accretion rate and with the X- ray emission 

1 Heinz fc Sunvaev 20031 iMerloni et al. 20031) as discussed 
for GX 339-4 i|Zdzia,rski et al. 20041) . However, detection 
of radio emission was reported du ring the TDS, for 
example in the BH XTE J1650-500 JCorbel et al. 20041) 
and the scenario increases in complexity (see Fender 
et al. 2005 for a complete review). More recently, Markoff 
et al. (2004) have proposed jet models where the synchro- 
self Compton or the external Comptonization radiation 
are the dominant processes generating X-ray spectra in 
BH binaries. These models seem to fit the experimental 
data as well as the thermal Comptonization models do. 

While we did not observe any significant high-energy 
tail during the PV-Phase observations, it appears that, in 
the SPI spectra of 2003 June, when the source was in the 
IS, some data points above 400 keV are not well described 
by the thermal Comptonization models (see Fig. U3 and EJ). 
All the tests we have performed using different SPI back- 
ground models and the Crab nebula spectrum indicate 
that this feature comes from the source. Similar excess 
over a Comptonization law has be en previously observed , 
in particular with CGRO/OSSE ^McConnell et al. 2002^1 
which detected a power law like component extending 
beyond 1 MeV during a (so classified) TDS of Cygnus X-l. 
In our case, the SPI data of the observed IS spectrum of 
Cygnus X-l imply a rather bright high-energy tail. The 
bolometric luminosity of the source, including this new 
component, is 1.65 times higher than the one obtained 
with the best-fit thermal current model. 

This kind of steep power law, without a high-energy 
break (at least up to 1 MeV) was modeled in the past 
including a non-thermal component in the accretion flow, 
with the so-called hybrid thermal/non-thermal models 
from Poutanen & Svensson (1996) and Coppi (1999). 
Using these same models, we obtain a better fit of our 
data, which indicate the necessity to include a non- 
thermal distribution of electrons, with a power law energy 
index between 2.4-3.7, depending on the model used. T p 
between 2 and 3 are expected from shock acceleration 
models. We find that most of the crucial parameters such 
as r, Q/2ir, Ep c centroid, EW or kT c are compatible with 
previous results reported by McConnell et al. (2002), 
Frontera et al. (2001) and Gierlihski et al. (1999). With 
the non-thermal EQPAIR model, we found an unabsorbed 
bolometric luminosity of 3.3 x 10 37 ergs s _1 (higher 
than the one obtained with the current thermal model). 
This value is ~ 1.9 times lower than the one observed by 
McConnell et al. (2002) in the (so classified) TDS: this, 
along with the differences seen in the lh/l s , ^nth/^h and T p 
values (see Section 3.2), definitely shows that our epoch 

2 observations did not happen during a TDS but rather 
in the IS, with an intermediate luminosity. 

Alternatively, Comptonization on a population of 
(thermalized) electrons with bulk motion, e.g., Titarchuk 
et al. (1997), Laurent & Titarchuk (1999), is sometimes 
invoked to explain the power law high-energy component. 
This could describe the observed high-energy emission 
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seen in epoch 2, as it predicts T > 2 or even softer, even 
if the lower energy cutoff should be around 100 keV. In 
addition to hybrid thermal/non-thermal models presented 
above, a stratified Comptonization region, providing a 
larger range of both electron temperatures and optical 
depths, could model the spectrum. Ling et al. (1997) 
reached the same conclusions based on Monte Carlo 
modeling of CGRO /BATSE spectra combined with (non 
contemporaneous) CGtfO/COMPTEL data (McConnell 
et al. 1994). Thermal gradients are incorporated into 
several other models, which then lead to the generation 
of a high-energy tail, e.g., Skibo & Dermer (1995), 
Chakrabarti & Titarchuk (1995) and Misra & Melia 
(1996). 



5. Summary and conclusions 

Using the broad-band capability of INTEGRAL, it has 
been possible to accumulate a large amount of simulta- 
neous data on Cygnus X-l between 5 keV-1 MeV and 
to follow its spectral evolution from 2002 November to 
2004 November. These data were helpful to characterize 
the evolution of the Comptonization parameters of the 
source correlated to the presence of a variable disc emis- 
sion, indicating transitions between the HS and a softer 
(Intermediate) state. We also observed the presence in 
the SPI data of a high-energy tail during the IS (or the 
FST), emerging from the Comptonization component be- 
tween 400 keV-1 MeV and probably associated with a non- 
thermal component. The extent to which the spectrum 
hardens at energies approaching 1 MeV has now become 
an important issue for theoretical modeling of the accre- 
tion processes and radiation mechanisms in BH binaries. 
Data from both IBIS and SPI instruments offer the best 
opportunity to define more precisely the high-energy X- 
ray binary spectra. We hope to further investigate this 
using INTEGRAL data from this source and other bright 
BH X-ray binaries. 
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